Extended MSSM Neutralinos as the Source of the PAMELA Positron Excess 
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We consider a scenario within the Minimal Supersymmetric Standard Model extended by a singlet 
chiral superfield, in which neutralino dark matter annihilates to light singlet-like Higgs bosons, which 
proceed to decay to either electron-positron or muon-antimuon pairs. Unlike neutralino annihilations 
in the MSSM, this model can provide a good fit to the PAMELA cosmic ray positron fraction excess. 
Furthermore, the singlet-like scalar Higgs can induce a large Sommerfeld enhancement and provide 
an annihilation rate sufficient to accommodate the observed positron excess. 
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The satellite-based experiment PAMELA has reported 
a cosmic ray positron fraction (defined as the ratio of 
positrons to electrons-plus-positrons) which rises rapidly 
between 10 GeV and 100 GeV [lj. This is in stark 
contract to the behavior predicted for positrons pro- 
duced through interactions of cosmic ray protons with 
the interstellar medium [2]. Although the origin of 
this positron excess is currently unknown, a number 
of plausible sources have been proposed, including pul- 
sars [3| , the acceleration ofpositron secondaries in cosmic 
ray acceleration regions [4, and dark matter annihila- 
tions [1, i, 0, !, ! Ell EElli or decays [H | . 



Efforts to explain these observations with annihilating 
dark matter face some challenges. In particular, the rapid 
rise of the PAMELA positron fraction appears to require 
a very hard injected spectrum, which in turn requires 
the responsible dark matter particle to annihilate primar- 
ily to charged leptons [5|, |g, [j], l8[ (for an exception, see 
Ref. |l5|). Furthermore, dark matter candidates which 
annihilate largely to quarks or gauge bosons are also pre- 
dicted to overproduce cosmic ray anti pro tons (ill . fl3 |. 



gamma rays, and synchrotron emission [13] if the overall 
annihilation rate is normalized to produce the positron 
fraction reported by PAMELA. Within the context of the 
Minimal Supersymmetric Standard Model (MSSM) , neu- 
tralino dark matter annihilates largely to final states con- 
sisting of heavy quarks or gauge and/or Higgs bosons 
As a result, such annihilations produce a relatively soft 
spectrum of cosmic ray positrons (l7j and are unable to 
provide a viable explanation for the PAMELA excess. 



The same conclusion is not necessarily reached in su- 
persymmetric models with an extended Higgs sector. Ex- 
tensions of the MSSM by a singlet chiral superfield are 
motivated in order to explain the size of the /i term [l8| , 
to raise the mass of the lightest CP even Higgs boson 
above the LEP II bound [19| , to reduce electroweak fine- 
tuning |20] , or to catalyze electroweak baryogenesis 21 1 . 



Such extensions are described by superpotential, 
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and soft Lagrangian, 
\m 2 s \S\ 2 + B S S 2 d 



\A x SH u H d + kA k S 3 + B.c. (2) 



Specific implementations of the singlet typically involve 
a subset of these terms. For example, the Next-to-MSSM 
(NMSSM) [IS] invokes a Z3 symmetry to remove all 
but the terms involving A, k, A\, and A K , whereas the 
Fat Higgs models [22[ have a dynamically generated su- 
perpotential utilizing the v 2 , term to drive electroweak 
symmetry-breaking, even when supersymmetry is unbro- 
ken. In this work, we do not wed ourselves to any one of 
these specific realizations, but find that our conclusions 
can hold for regions of parameter space in any of them. 

The additional singlet results in an extra neutralino, 
and two scalars (one CP even and the other odd) in the 
spectrum. The parameter space of interest here are re- 
gions in which the lightest neutralino is largely singlino 
and the lightest CP even and odd scalars are largely 
singlets. In this region, the singlet constitutes a kind 



of hidden sector [23( which mixes with the MSSM Hig- 
gses through electroweak symmetry-breaking. A mostly 
singlino neutralino is relatively simple to arrange in the 
limit in which (fj, s + K(S)) < M u M 2 , (n + X{S)). Alight 
and mostly singlet pseudoscalar arises naturally when 
A\, A K , m|, Bs, v 2 , and us are small, because in this 
limit it is the the pseudo-goldstone boson of an explicitly 
broken U(l) symmetry. As we will see below, it is also 
preferable to have a light and mostly singlet scalar boson. 
Unlike the singlino and pseudoscalar, this requires some 
engineering of parameters, but the tunings involved are 
relatively modest, at the 10 — 20% level. 

In the limit of a singlino-like lightest neutralino to- 
gether with light singlet-like scalar (h) and pseudoscalar 
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(a) Higgs bosons, neutralino annihilations proceeds dom- 
inantly to a ah final state through i/u-channel neutralino 
exchange and s-channel a exchange, with a low velocity 
cross section given (in the limit of m x o 3> m ai rrih) by: 
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where I), xx and T axx are the Higgs couplings to the neu- 
tralino, 3h aa is the coupling between the Higgs bosons, 
and v is the relative velocity between the WIMPs. In the 
limit of singlet-like Higgs bosons and a singlino-like neu- 



tralino, these couplings reduce to T axx m Tj 



-V2k, 

and ghaa ~ k(3A k — /is — k(S)). Singlinos can also anni- 
hilate to hh or aa final states, but the rates for these pro- 
cesses are suppressed by v 2 . In the regime ghaa/mx ~ ^ 
the cross section depends only on n and m x o , and results 
in a thermal relic abundance of approximately, 



200 GeV 



(4) 



If rrih > 2m a , the scalar Higgs bosons produced in 
the annihilations will decay dominantly to a pair of the 
pseduoscalar Higgs bosons, leading to a 3a final state. 
The pseudoscalar Higgs bosons decay via a small mix- 
ing angle with MSSM Higgs bosons, and may be some- 
what long lived. Typically, these decays proceed to 
the heaviest kinematically available fermions. As we 
are interested in scenarios in which neutralino annihi- 
lations yield mostly charged leptons, we focus on the 
case m a < GeV, which leads to the production of ei- 
ther muon pairs (2m ^ < m a < 2m 7r ) or electron-positron 
pairs (2m M > m a ). 

To calculate the spectrum of positrons and electrons 
in the cosmic ray spectrum, we model the diffusion and 
energy losses of such particles. This is done by solving 
the steady-state propagation equation 24j |: 
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where dn e /dE e is the number density of elec- 
trons/positrons per unit energy, K(E e ) is the diffusion 
coefficient, and b(E e ,x) is the energy loss rate. The 
source term, Q(E e ,x), reflects the mass, annihilation 
cross section, dominant annihilation modes, and distri- 
bution of dark matter in the Galaxy. We adopt a dif- 
fusion coefficient of K(E e ) = 6.04 x 10 28 (£' e /4 GeV) 0,41 
cm 2 /s, and boundary conditions corresponding to a disk 
of 5 kiloparsecs half-thickness, which yields the best-fit to 
the current body of cosmic ray data (stable and unstable 
primary-to-secondary nuclei ratios) [lH and an energy 
loss rate of b(E e ) = 10~ 16 (i? e /l GeV) 2 GeV/s, resulting 
from inverse Compton and synchrotron processes. For 
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FIG. 1: The cosmic ray positron fraction resulting from neu- 
tralino annihilations in several selected scenarios. In each 
case, we consider the channel x°X° ~* followed by h — * aa, 
and a —> e + e~ (top) or fi + fi~ (bottom). In the upper (lower) 
frame, we have used mi = 3 GeV (10 GeV), although the 
precise value of this mass has only a small effect of the shape 
of the electron/positron spectrum. In each case, we have nor- 
malized the annihilation rate to accommodate the PAMELA 
data. The dot-dashed line denotes the prediction from astro- 
physical secondary production alone. 



the source term, we adopt a Navarro-Frenk- White halo 
profile to describe the dark matter distribution in the 
Milky Way. To convert the positron/electron spectrum 
to a positron fraction, we use the primary and secondary 
cosmic ray spectra as described in Ref . |2j . 

In Fig. [TJ we show the cosmic ray positron fraction re- 
sulting from neutralino annihilations for several choices of 
the neutralino and pseudoscalar Higgs masses. In the top 
frame, we consider Higgs decays that produce electron- 
positron pairs (m a < 2m M ), whereas in the lower frame, 
the decays proceed to muon pairs (m a > 2m fl ). The 
mass of the scalar Higgs only mildly impacts the result- 
ing spectrum of positrons. 

In each case shown, we have normalized the neutralino 
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Required BF 


Corresponding 


130 GeV 


200 MeV 


40 


< 2 GeV 


200 GeV 


1.22 MeV 


11 


~ 1-4 GeV 


200 GeV 


200 MeV 


45 


< 4 GeV 


600 GeV 


212 MeV 


90 


<, 30 GeV 


600 GeV 


250 MeV 


100 


<> 30 GeV 



TABLE I: The boost factor to the neutralino annihilation rate 
required to produce the PAMELA positron excess, for sev- 
eral values of the neutralino and singlet pseudoscalar Higgs 
masses, and the approximate value (or range) of the singlet 
scalar Higgs mass required to generate the required boost fac- 
tor through the Sommerfeld effect. 



annihilation rate to accommodate the PAMELA signal. 
This rate scales with the neutralino annihilation cross 
section and with the square of the local dark matter den- 
sity. The annihilation rate may also be enhanced as a 
result of inhomogeneities in the dark matter distribution 
(ie. clumps or other substructures). Using a local dark 
matter density of 0.4GeV/cm 3 and an annihilation cross 
section of av = 3 x 1CP 26 cm 3 /s (required to thermally 
generate the observed abundance of dark matter) , we find 
that the annihilation rate must be enhanced by a factor 
of ^10-50 for annihilations to electrons or ^100 for an- 
nihilations to muons. 

Fortunately, the rate of neutralino annihilations in the 
Galactic Halo can be naturally enhanced in this scenario 
through the Sommerfeld effect generated by the light 
scalar Higgs boson. Although the resonance structure 
of the Sommerfeld effect can be complex, very roughly 
speaking, if the light Higgs scalar is lighter than nih 5a 
K 2 m x o /4-7T, then for a velocity dispersion expected for the 
Galactic Halo (a ss 150 km/s), the annihilation rate will 
be enhanced by a factor of about S ~ n 2 m x o/4irmh- 
More details can be found i.e. in the appendices of 
Ref. 0. 

In Table I, we list the enhancement to the neutralino 
annihilation rate (relative to that obtained for av = 
3x 10 -26 cm 3 /s, no significant substructure, and no Som- 
merfeld effect) required to normalize the positron fraction 
to the PAMELA observations. For each parameter set, 
we determine (assuming k is such that the correct ther- 
mal relic density is predicted, as per Eq. [4J the approx- 
imate value of the scalar Higgs mass that would lead to 
a Sommerfeld enhancement providing the required boost 
factor to explain the PAMELA data. We find that for 
m x o ~ 200 GeV (600 GeV), we must require rrih 5o 4 GeV 
(< 30 GeV) in order to generate the desired annihilation 
rate. 

Light scalars and pseudoscalars must be mostly sin- 
glets in order to have escaped detection by previous ex- 
periments. A 30 GeV scalar should be more than about 
90% singlet, or it would have been discovered by LEP 
II in its Standard Model Higgs search 27j. Tevatron 



data provides an important constraint through produc- 
tion of a scalar Higgs which decays through two light 
pseudoscalars into a four muon final state [28|. The 
null result of this search will be satisfied for the lightest 
CP even Higgs provided it has already escaped the LEP 
bound, and further requires the heavier CP even MSSM 
Higgs bosons to have masses greater than about 135 GeV 
to suppress the decay into pseudoscalars. For the very 
light a masses considered here, there is also the possibility 
of observing rare decays such as T(3s) — > ja 
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30j and K+ -> ir+a -> ir+fi+fj,- [3l|,l32j. Both mea- 
surements are consistent with an a which is at least 90% 
singlet, provided tan/3 is of order one. 

Thus far, we have not addressed the electron (plus 
positron) spectrum as measured by ATIC [33| . and 
more recently by the Fermi Gamma Ray Space Tele- 
scope (FGST) [3J|. The sharp, edge- like feature at 
~600 GeV reported by ATIC could easily be accommo- 
dated in the scenario discussed here for m x ~ 600 GeV, 
m a < 2m M , and < 30 GeV. The spectrum newly 
reported by FGST could also potentially be accommo- 
dated, but would require a multi-TeV mass for the light- 
est neutralino [3q . 

In summary, we have presented a scenario within the 
context of the Minimal Supersymmetric Standard Model 
extended by a gauge singlet in which the lightest neu- 
tralino annihilates to light singlet- like Higgs bosons which 
proceed to decay to either electron-positron or muon- 
antimuon pairs, leading to a cosmic ray positron frac- 
tion consistent with observations of the PAMELA exper- 
iment. Furthermore, the annihilation rate of the neu- 
tralinos in the Galactic Halo can be strongly enhanced 
by the Sommerfeld effect in this model. No astrophysical 
boost factors are required to obtain the positron fraction 
observed by PAMELA. 
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Note: As this letter was being finalized, a related study 
appeared Xshi which discusses NMSSM neutralinos as a 
source of the PAMELA excess, but in a considerably dif- 
ferent region of parameter space. 
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